Decoherence of localized spins interacting via RKKY interaction 
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We theoretically study decoherence of two localized spins interacting via the RKKY interaction in 
one-, two-, and three-dimensional electron gas. We derive the kinetic equation for the reduced density 
matrix of the localized spins and show that energy relaxation caused by singlet-triplet transition is 
suppressed when the RKKY interaction is ferromagnetic. We also estimate the decoherence time of 
the system consisting of two quantum dots embedded in a two dimensional electron gas. 
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Quantum computation and quantum information are 
emerging research fields of physics, technology and infor- 
mation sciences^. The elementary units in most quan- 
tum computation and quantum information schemes are 
a quantum bit and a quantum gate. Because of its 
scalability and relatively long coherence timeQ], solid 
state device with localized spins is considered as one of 
the promising candidates for these quantum devices^- 
KaneQ proposed the system of nuclear spins of phos- 
phorus donors in a silicon heterostructure with direct 
exchange interaction between electrons localized at the 
donors. Mozyrsky et al. 0, |(| proposed the qubits of nu- 
clear spins in quantum-Hall system in which the gate 
operation is realized by the indirect exchange interac- 
tion via virtually excited spin waves. Recently Craig 
et a/.@, 0] observed the RKKY coupling of semi- 
conductor quantum dots, which proved that the RKKY 
interaction can be used as a quantum gate consist- 
ing of localized spins in semiconductor quantum dots. 
Several theoretical papers studying the RKKY interac- 
tion in such semiconductor nanostructures have been 
published[l(J El 13 El 13 E3 

One of the major obstacle to realizing quantum com- 
putation is decoherence EE> El- Since the RKKY inter- 
action is mediated by the electron gas, the particle-hole 
excitations act as an environment (Fermion bath) [TM [Tgt 
I20L l2l| . It is therefore important to clarify the effects of 
the Fermion bath on the dynamics of the qubits consist- 
ing of the localized spins. 

In this paper, we study the dynamics of the localized 
spins interacting via the RKKY interaction by using the 
kinetic equation of the reduced density matrix. We find 
the term intrinsic to the RKKY interaction appears in the 
kinetic equation, which is an oscillating function of the 
distance between the localized spins with the same pe- 
riod as the RKKY interaction. We show that the energy 
relaxation due to the singlet-triplet transition is strongly 
suppressed by this term. We also discuss the physical re- 
alization of a quantum gate using the RKKY interaction 
and estimate the decoherence time of the system consist- 
ing of two quantum dots embedded in a two dimensional 
electron gas(2DEG). 

We consider the system consisting of two localized 
spins embedded in a one-, two-, or three-dimensional elec- 
tron gas. The Hamiltonian, H — Hs + H c + Hi nt , com- 



prises the localized spin part, Hs, the conduction electron 
part, H c , and the s-d interaction, H mi . We assume that 
there is no external magnetic field and we set H$ = 0. In 
the second quantization representation, H c and Hi nt are 
expressed as 

H c = ^2e k c{ a c kry , (1) 

ka 
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where ct and Cka are creation and annihilation opera- 
tors of an electron with wavenumber vector k and spin 
a. Here Sk is an energy of conduction electrons, J is a 
coupling constant of s-d interaction, and R p (p — 1,2) 
represents the position of the localized spin S p . We as- 
sume that the coupling constant J is so small that we 
can treat -Hint as a perturbation on H c . 

The dynamics of the two localized spins is described 
by the reduced density matrix p{t) = tr c {ptot(£)}, where 
Ptot is the density matrix for the total system and tr c 
means trace over the degrees of freedom of the conduction 
electrons. The reduced density matrix p(t) obeys the 
following kinetic equation |22l 123. l24l |25| : 

j fP (t) = -- [H s , p(t)] + J* dt'E(f - t')p(t'). (3) 

The self-energy E(t — t') is a super-operator acting on 
pit). We assume that at the initial time t = the system 
of the localized spins and the system of the conduction 
electrons are decoupled, and the conduction electrons are 
in thermal equilibrium: ptot(0) = p{Q)®p^ q , where p^ q = 

e -(H a -nN)/k B T / tlce -(H a -nN)/k B T_ 

We carry out a second-order perturbative calculation 
for the self-energy E(t — t') in Eq.©. We assume that 
the reduced density matrix p{t) varies slowly compared 
to the lifetime of particle-hole excitations, and we make 
the Markov approximation in Eq. |J3J| . A straight forward 
calculation gives 

^-p(t) = ( Hrkky + V p + v «* J ( 4 ) 
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TABLE I: The coefficients a d and /3d, and the range functions Fd(x) and Gd(x) for d-dimensional electron gas (d = 1,2,3). 
si(x) = — dtsmt/t is the sine integral function. J n is the Bessel function and N n is the Neumann function of order n. 
zo = 2.40 is the first zero-point of Jo. 



where Wrkky, 2?p (p — 1,2), and 2? ex are the super- 
operators defined as follows. 

The first term of r.h.s. in Eq.Q is defined as 

HrkkypM = -z^p [ Sl . s 2l p (t)] , (5) 

which represents the coherent time evolution of the two 
localized spins interacting via the RKKY interaction^ 
Hi! HH < The effective coupling constant Jrkky, for 
d-dimensional system (d = 1,2,3) is given by Jrkky = 
adPdEpF d {2k F R), where R = \Ri — R 2 \ is the distance 
between two localized spins, Ep is the Fermi energy, and 
kp is the Fermi wavenumber. Here ad, (3d, and the range 
function Fd(x) are given in Table [I] The range functions 
Fd(2kpR) for d = 1, 2, and 3 are plotted by dotted lines 
in Figs^a), (b), and (c), respectively The RKKY inter- 
action is produced by virtually excited quantum states 
within a wide energy window from the bottom of energy 
band to the Fermi surface. Therefore, the strength of the 
RKKY interaction, Jrkky, is proportional to the Fermi 
energy Ep and is oscillating function of R with the period 
of half of the Fermi wavelength, w/kp. 

The second term of r.h.s. in Eq.QJ is defined as 

£ v pP {t) = -{ £ (|p(*) - (s+ P (t)s- 

p=l,2 p=l,2 (6) 

+ s-p(t)s+ + 2S z pP (t)s;)), 

where the coefficient 7 is expressed as 7 = a^fcaT This 
term describes the usual decoherence of a localized spin 
S p interacting with the Fcrmion bath, which is known 
as the Korringa relaxation |2^|. Since thermally excited 
particle-hole pairs cause the Korringa relaxation, 7 is 
proportional to the temperature T. One can easily show 
that this term causes energy relaxation due to transition 
between singlet and triplet states and dephasing due to 
transition among three degenerate triplet states. 

The last term of r.h.s. in Eq.Q, which is intrinsic to 



the system with the RKKY interaction, is defined as 

v e * P {t) = ■ s 2 , P (t)} Us+ P (t)s 2 - 

h V 2 (7) 

+ S^p(t)S+ + 2S z lP {t)S z 2 + (spin 1 «-> 2))J. 

This term exists only when the two localized spins inter- 
act with each other via conduction electrons. The coef- 
ficient 7 CX is expressed as 7 ex = otdkBTGd{2kpR). The 
Range function Gd(2kpR) is given in tabled and plotted 
by solid lines in Figsnja),(b), and (c). As we shall show 
later, this term suppresses the energy relaxation caused 
by singlet-triplet transitions due to the term T> p . The ori- 
gin of this term is interference among thermally exited 
particle-hole pairs near the Fermi surface. Therefore, the 
range function Gd(2kpR) oscillates with the same period 
as Fd(2kpR). Reflecting the difference in physical mech- 
anisms behind them, the amplitude of Gd{2kpR) decays 
as 1/R d ^ 1 while that of Fd(2kpR) decays as slowly as 
1/R d . Especially for the one-dimensional electron gas, 
the amplitude of the oscillation of Gd(2kpR) does not 
decay as R increases. 

Next, in order to capture the physical meaning of the 
term V cx p(t), we examine the time evolution of the ex- 
pectation value (Si ■ S 2 ) = tr{p(t)Si ■ S 2 }, which is pro- 
portional to the energy expectation value of the RKKY 
interaction. By using Eqs.Q}-©, we can show that it 
decays exponentially as (Si ■ S 2 ) cx exp(— 4(7 — 7 cx )t//i). 
Then the characteristic time of energy relaxation of the 
RKKY interaction is given by trkky = ^/4(7 — 7cx)- 
Since the term V cx p(t) is proportional to 7 0X , we can say 
that the term V ex p(t) tends to suppress energy relax- 
ation. Note that 7 ex does not exceeds 7. Especially when 
7 0x = 7, the singlet-triplet transitions due to the term 
T> p does not occur and energy relaxation of the RKKY 
interaction is completely suppressed. However, the tran- 
sitions among the triplet states are not suppressed and 
dephasing rate remains finite even if 7 CX = 7. As shown 
in FigsQJa), (b), and (c), the range function Gd(2kpR) 
takes its maximal values when the Fd(2kpR) is negative. 
Consequently the energy relaxation of the RKKY inter- 
action is suppressed when the RKKY interaction forms 
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FIG. 1: The range functions for one-dimensional system 
Fi(2kpR) and Gi(2kpR) are plotted as a function of 2kpR 
in Panel (a). Panels (b) and (c) are the same plot for two- 
and three-dimensional systems, respectively. 



ferromagnetic coupling (Jrkky < 0). 

Figure Ola) shows the time evolution of the diagonal 
element Pn,n(t) of the reduced density matrix. The ini- 
tial state is taken to be |f|) and the time scale is nor- 
malized by hj Jrkky- The diagonal element Ptx.TlW 
represents the occupation probability of the state \H). 
Without decoherence, 7 = 7 GX = 0, the quantum state 
of the two-spin system oscillates between |fj.) and |J.|) 
coherently and p-fj,,|j,(i) shows a clear oscillation like a 
trigonometric function as indicated by the thin solid line 
in Fig|2fa). In the case of 7 = 0.05 Jrkky and 7 CX = 0, 
where each spins decohere independently /0|x,ti W sriOWS 
damped oscillation due to V p (thick solid line). The 
dotted, dashed, and dot-dashed lines in Fig. EJa) show 
the results for 7 CX =0.025, 0.045, and 0.05Jrkky, respec- 
tively. The value of 7 ex can be controlled by changing the 
distance between localized spins, R. One can see that the 
oscillation amplitude of p-\i,n{t) decays more slowly as 
7 ex increases since the energy relaxation of the RKKY 
interaction is suppressed by 7 CX . Although energy relax- 
ation of the RKKY interaction is completely suppressed 
when 7 0X — 7, dephasing due to the transition among the 
triplet states occurs as indicated by the dot-dashed line 
in Fig. Ha). 

In Fig. |2Ib) we plot the off-diagonal element Pi 141 (t) 
with the initial state (\H) — \[]))/^/2 (the singlet state, 
i.e., the maximally entangled state of the localized spins). 
Because the singlet state is one of the energy eigenstates 
of the RKKY interaction, Pn,n{t) is the conserved quan- 
tity when 7 = 7 CX = 0. If 7 takes a finite value, the sin- 
glet state changes to the mixed state and Pu,ll(t) decays 
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FIG. 2: (a) The time evolution of Ptl.uW- The initial 
state is taken to be and the time is normalized by 

h/JnKKY- (b) The same plot for PTI.nM with the initial 
state (HI) — \l1})/\/2. For both panels, thin solid line rep- 
resents the result for 7 = 7 CX = while all the other lines 
are for 7 = 0.05 Jrkky- The parameter 7 CX is taken to be 
0, 0.025, 0.045, and 0.05Jrkky for the thick solid, dotted, 
dashed, dot-dashed lines, respectively. 



exponentially to zero. As shown in Fig|3 (b), PtL-ltW 
decays more slowly as 7 GX increases, which means that 
7 CX suppresses decoherence of the singlet state. Espe- 
cially in the case of 7 = 7 GX , the singlet state remains as 
it is forever and the off-diagonal element Pti,iT (^) never 
decays although the localized spins always interact with 
the Fermion bath of conduction electrons. One can eas- 
ily show that for the singlet state the contributions of 
T3 C xp{t) and T> p p{t) in Eq.Q) cancel out when 7 = 7 CX . 

Finally, we would like to discuss the physical realiza- 
tion of a quantum gate using the RKKY interaction and 
estimate the decoherence time. Let us consider the sys- 
tem consisting of two quantum dots embedded in a 2DEG 
of GaAs/AlGaAs heterostructure. We assume that the 
electron density of the 2DEG is n 2D = 7.3 x 10 n cm~ 2 , 
the Fermi energy and the Fermi wavelength are Ep ~ 
26meV and 2ir jkp ~ 29.4nm. We also assume that the 
charging energy of the quantum dot is U = 1.9meV, the 
Lorentzian broadening of the localized-state energy with 
full width at half maximum (FWHM) is T = 295/xeV, 
and the temperature is T = lOOmK [H H3- Each 
quantum dot contains a single localized spin which acts 
as a qubit. The coupling between localized spins, i.e., 
the quantum gate operation is controlled by applying a 
gate voltage to 2DEG. For example the VSWAP gate 
operation^, l32| . which is known as a universal 2-bit 
quantum gate operation, can be carried out if strength 
of the RKKY interaction Jrkky is controlled such that 
J dt Jrkky (t)/h = tt/2. We suppose that the inter dot 
distance is R = z^/Hf ~ llnm (S2 — 2.40 is the first zero- 
point of the Bessel function Jo), that is, 7 CX = and the 
energy relaxation is not suppressed. The strength of the 
RKKY interaction takes Jrkky ~ 3.4/ieV and the oper- 
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ation time defined as r op = 7t/i/2Jrkky is ~ 0.3ns. Since 
7 0X is assumed to be zero, the decoherence time is deter- 
mined only by 7 and is estimated as Td cc = V7 ~ 62ns. 
Therefore, about 200 times coherent VSWAP operation 
can be achieved. 

The decoherence of a localized spin due to phonon scat- 
tering is suppressed in a small quantum dot and long 
coherence time (order of millisecond) of a single spin 
was observed in self-assembled semiconductor quantum 
dots0. Therefore, the decoherence due to the Fermion 
bath of conduction electrons is dominant in this system. 

More coherent quantum gate operation is available 
if we locate the quantum dots very close to each 
other(kpR <C 1). In this case the value of 7 CX ap- 
proaches 7 and therefore we can take an advantage of 
the suppression of the energy relaxation. Furthermore 
we can also obtain very strong ferromagnetic RKKY 
interactional]] ■ Experimentally quantum dot arrays 
with 3-10nm diameter quantum dots have already been 
fabricated |33l I34L l35| . We believe that the system we con- 



sider can be realized by using such quantum dot arrays 
with a few nm inter dot distance. 

In conclusion, we have derived the kinetic equation for 
the system of two localized spins embedded in an electron 
gas and have shown that particlc-hole excitations in the 
electron gas cause not only the RKKY interaction but 
also decoherence of the two-spin system. We also show 
that the strength of decoherence as well as the RKKY in- 
teraction strongly depends on the distance between two 
spins, and energy relaxation due to singlet-triplet transi- 
tion is suppressed when the RKKY interaction is ferro- 
magnetic (Jrkky < 0). We also estimate the decoher- 
ence time of the system consisting of two quantum dots 
embedded in a two dimensional electron gas(2DEG) to 
be Tdcc ~ 60ns within which ~ 200 coherent VSWAP 
operations can be achieved. 
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